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I. INTRODUCTION 

There is a variety of physical phenomena in which the 
processes of diffusion, coagulation, adsorption and des- 
orption play an important role. For example, submono- 
layer epitaxial thin film growth involves deposition of 
atoms onto a substrate and diffusion of these atoms lead- 
ing to their aggregation into islands of increasing size . 
A second example is river networks which have been mod- 
eled by aggregating masses in a steady influx of particles 
0, El ■ Further examples include aerosols and clouds , 
colloids and polymerization [||. 

A simple lattice model that incorporates the above pro- 
cesses is the In-Out model 0,111 in which diffusing point 
size particles on a lattice coagulate together on contact 
forming particles of larger mass. In addition unit mass 
is input uniformly at rate q while unit mass evaporates 
from an existing mass at rate p. The competition be- 
tween adsorption and desorption results in a nonequi- 
librium phase transition between a phase in which the 
average mass in the system is finite to one in which it 
increases linearly with time. A quantity that captures 
the features of the steady state is the mass distribution 
P(m,t). P(m,t) is the probability that a randomly cho- 
sen site has mass m at time t. For fixed desorption rate 
p, the distribution P(m, t) for large times was shown to 
change from an exponential distribution at small values 
of q to a power law m~ Tc at q = q c (p) and to a differ- 
ent power law m~ T for q > q c (p). Near the transition 
point, P(m, q — q c , t) was seen to have the scaling form 
P(m, q — q c , t) ~ ni~ Tc Y [m(q — q c )^, mt~ a ~\ . The values 
of these exponents in high dimensions were calculated us- 
ing a mean field approximation. In one dimension, they 
were determined using monte carlo simulations. The nu- 
merical values in one dimension were significantly dif- 
ferent from the mean field results. In this paper, these 
exponents are calculated in all dimensions. The model 
is also extended to one in which particles diffuse with a 
mass dependent rate with /i > 0. The critical ex- 
ponents for this more general model are also calculated. 

Related models have been studied in the context of 
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nonequilibrium wetting [{j, ll(| . If the mass m in the In- 
Out model is identified with the height of a substrate, 
then the transition observed is qualitatively similar to 
nonequilibrium wetting transitions. In these models, the 
system undergoes a transition from a phase in which the 
interface is smooth to one in which it is rough. The expo- 
nents describing these transitions have been found to be 
related to some underlying contact process undergoing an 
absorbing to active transition. The In-Out model stud- 
ied in this paper differs from these models by the lack 
of a surface tension term which tries to smoothen out 
the interface, and thus belongs to a different universality 
class. 

The rest of the paper is organized as follows. In SeclTTl 
the model is defined, the results obtained in Ref. 8] for 
the \i = case are reviewed, and the results of this paper 
are summarized. In Sec. IIIII a scaling relation is derived 
between the critical exponents. In Sec. If VI the phase in 
which the mean mass increases linearly with time is fully 
characterized. In Sec. |3 exactly solvable limits of the 
model, namely /i = I case, zero dimensions and mean 
field solution, are discussed. In Sec. IVII the exponents 
for arbitrary /i are derived in all dimensions with the 
help of an assumption. The results are compared with 
monte carlo simulations in one dimension. In Sec. IVIII 
the results of the model are compared with results of 
related models. Section IVIIII contains a summary and 
concluding remarks. The appendices contain the details 
of the calculations. 



II. MODEL AND RESULTS 

A. Definition 

For simplicity, we define the model on a one dimen- 
sional lattice with periodic boundary conditions; gener- 
alizations to higher dimensions is straight forward. Each 
site i of the lattice has a non-negative integer mass vari- 
able rrii > 0. Given a certain configuration of masses 
at time t, the system evolves in an infinitesimal time dt 
as follows. A site i is chosen at random (with proba- 
bility dt), and then the following events can occur, (i) 
Adsorption: with probability q/(p + q + 1), unit mass is 
adsorbed at site i; thus m, — * rrii + 1. (ii) Desorption: 
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if the mass rrii is greater than zero, then with proba- 
bility p/(p + q + 1), unit mass is desorbed from site i; 
thus m,i — > — 1 provided m, > 0. (iii) Diffusion and 
aggregation: if the mass m t is greater than zero, then 
with probability m i ^ / (p + q + 1) the mass rrii moves to 
a randomly chosen nearest neighbor. If the target site 
already happens to have some mass, then the total mass 
just adds up; thus mi — * and mm — * mj±i + nij. The 
initial condition is chosen to be to be one in which all 
sites have mass zero. The model has three parameters, 
p,q,fi. 



B. Review of results for /j, — 

For ii = 0, when all particles diffuse at the same rate, 
the single site mass distribution at time t, P(m, t) was 
determined in large dimensions using a mean field ap- 
proximation and in one dimension using monte carlo 
simulations @. It was shown that when the adsorp- 
tion rate q was increased keeping the desorption rate p 
fixed, the system undergoes a nonequilibrium phase tran- 
sition across a critical line q c {p) from a phase in which 
P(m) = limj^oo P(m, t) has an exponential tail to one 
in which it has an algebraic tail for large mass; i.e, 




when q < q c , 
when q = q c , 
when q > q c , 



(1) 



where m* is a q dependent cutoff, and r and r c are ex- 
ponents characterizing the power law decay. In addition, 
it was argued that as a function of the small deviation 
q = q — q c , and large time t, P(m, q, t) displays the scaling 
form 



P(m,q,t)^J-Y(mq' t ',^ 
m T " \ t a 



(2) 



in terms of three unknown exponents 4>, a, t c , and the 
two variable scaling function Y. The the three phases will 
be called as the exponential phase (q < q c ), the critical 
phase (q = q c ) and the growing phase (q > q c ). 

Of interest are two more exponents. The mass at each 
site could be thought of as representing the height of an 
interface. In that case, two quantities of interest are the 
average velocity of the interface F = d(m)/dt and the 
fluctuations of the interface about its mean which would 
be dominated by (m 2 ), where (. . .) denotes a spatial av- 
erage. In the exponential phase, (m) is finite and hence 
F = 0. In the growing phase, the interface has a finite 
velocity, and the velocity increases from zero as F ~ q e , 
where 9 is an exponent. At the critical point (m 2 ) ~ t@ . 
Using the scaling form Eq. it is straightforward to 
derive Hi 



9 = <f> [1 - a(2 - T C )] /a, 
= a(3 - t c ). 



(3) 
(4) 



The model when il = was studied using a mean field 
approximation that ignored the spatial correlations be- 
tween masses at different sites jg. It was shown that 
4> = 1, a = 2/3 and t c = 5/2 when fi = 0. Corre- 
spondingly, 9 = 2 and f3 = 1/3. In one dimension, the 
exponents when /i = were determined numerically to 
be t c w 1.83, a w 0.61, <j) w 1.01, 9 w 1.47, and (3 w 0.71. 



C. Summary of results 

In this paper, the In-Out model is studied for /i > 0. 
Using scaling arguments, a relation is derived between 
the exponents a and r c , thus reducing the number of 
unknown exponents from three to two. In particular, it 
is shown that 



a(jid + 2t c - 2) = d, d< 2. 
The exponent t c in d < 2 is shown to be 

= ) 2(d+2) M - u > 
c j(3_i^±2 0</i<2 . 

The exponent a in d < 2 is shown to be 



(5) 



(6) 



(7) 



The exponent <f> in d < 2 is calculated for = and 
M =l: 



1 M = 0, 

2 



M = l. 



(8) 



In the growing phase, the exponent t characterizing the 
power law decay of the mass distribution is shown to be 



(2 - n)d + 2 
d+2 



(9) 



In one dimension, when fj, = 0, the exponents reduce 
to t c = 11/6, a = 3/5 and 0=1. This is in vervgood 
agreement with the numerical results seen in Ref. [8j (see 
Sec. lIIB)l . In dimensions greater than two, the exponents 
take on their mean field value, obtained by setting d to 
2 in the above equations. 



III. SCALING RELATION BETWEEN a AND r c 

In this section, a relation between a and t c is obtained 
from scaling arguments for all /i. The dependence of the 
largest mass in the system M t on t can be obtained by 
the catchment area argument as follows. Due to diffusion, 
the mass M t would sweep out an area Lf in time t, where 
Lt is the typical length scale in the system. In addition 
to the mass contained in this area, M t also increases due 
to the average flux F = d(m)/dt. Thus, 



M t - L d t Ft. 



(10) 
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The typical length L t arises from diffusion: L t ~ (Dt) 1 / 2 . 
Substituting M t -/i for the diffusion rate D, and using 
F ~ t Q ( 2 ~ r c)- 1 from the scaling relation Eq. |J2J, we ob- 
tain M t ~ t(d+2a(2-r c ))/(2 +Al d)^ But; Dy definition (see 
Eq. J5J) M t ~ i Q . Equating the exponents, we obtain 

a(fid + 2t c - 2) = d. (11) 

Thus, the number of unknown exponents reduces from 
three to two. The above scaling arguments are valid only 
when L t increases as a positive power of t. This restric- 
tion translates to the condition afi < 1. 



implies that Y(x,y) ~ x Ta T fijj/x 1 a ^ s ) when x >■ 1. 
Thus 

f^-V-nfi!')' (15) 

where 7 is a crossover exponent. To make a comparison 
with Eq. (|14fl , one has to make the identification F ^ q e . 
Using Eqs. ©, l|ll|) and i|13|l . it is easy to show that 

~</>{t c -t) ^ pd/(d+2) _ Also, 

7 = 20/(<Z + 2). (16) 



IV. THE GROWING PHASE (q > g c ) 

In this section, the behavior of P(m) in the growing 
phase q > q c is discussed. The exponent r (as defined 
in Eq. is expected to be independent of the precise 
values of q and p as long as we are above the critical 
threshold y|. To obtain r, the convenient limit p = 
and q arbitrary may be studied. Different aspects of this 
limiting case have been studied in the context of river 
networks, self organized criticality and epitaxial growth 
IH 13 13 III El E3 We give a short derivation of the 
exponent r. In this limiting case of only adsorption, it is 
known that P(m, t) has the scaling form 



V. 



SOLVABLE LIMITS 



P(m, t) ~ m 



(12) 



where the scaling function f(x) tends to a constant (for 
fi < 1) for small values of x and decays exponentially for 
large values of x. Since there is a constant influx F of 
particles, (to) = Ft. Therefore 5(2 — r) = 1. To obtain a 
second relation between the exponents, note that Eq. i|ll|) 
is valid when r c is replaced by r and a by 6. Solving these 
two exponent equalities, one obtains 



In this section, we examine limiting cases of the model 
which are analytically tractable. For the sake of continu- 
ity of argument, the details of the calculation are deferred 
to the appendices. 



A. Solution for fj, — 1 

The special case when a mass m diffuses as to -1 can be 
solved by examining the time evolution of the two point 
correlations. When fj, = 1, certain simplifications occur. 
We refer to Appendix lAl for details. It is shown that the 
critical q c at which the mean mass increases with time is 



dp 2 g{p), 



(17) 



where 

g(p) 



2n dh 
„ 2^ 



271 



dkri 



2tt (l+p)d-£i=i cos ( fe i 



(18) 

The exponents for /i = 1 is shown to be (see Appendix Ia"|) 



S = 



(2- /i)d + 2 
d+2 ' 
d + 2 



fid + 2' 



(13) 



Equation (|13|l is valid when p < 1 and d < 2. For d > 2, 
the mean field results are correct. The fi — 0, results 
were obtained earlier in 

[3 13 G3 For n > 0, the 

one and two dimensional results were obtained earlier 
[3, 13|- The dependence of P(m, t) on the flux F can 
now be incorporated into Eq. (|12|) by simple dimensional 
arguments: 



F d + 2 



TO' 



(14) 



where r and 6 are as in Eq. i|13|) . 

The two variable scaling function Y(x, y) in Eq. (0) 
should be such that when x 1 (or mq^ 1), it reduces 
to the one variable scaling function / in Eq. (|14|) . This 



2d + 2 
d + 2 ' 
d + 2 



d 

d- 



(19) 
(20) 
(21) 



Solving for 9 from Eq. ©, we obtain 6 = {d + 2)/2. 



B. Mean field theory 

The exponents may be computed in large dimensions 
by a mean field analysis. This approximation involves 
ignoring the correlations between masses at two differ- 
ent sites. The details of the calculation are presented in 
Appendix The results are 



0, 



f < n < 2. 



(22) 
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The exponent a in the mean field equals 

M = 0, 



lh 0<M<2. 



(23) 



The exponent <f> could be computed only for ji = and 
M =l: 



rl m = 0, 
9 12 M=l- 

Using Eqs. and (|16l) . we obtain 
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2 /i 
1 M 



0,1, 
0,1. 



(24) 



(25) 
(26) 



Comparison with the exact solution for fj, = 1 or calcu- 
lating the dimension d when the mean field exponents 
satisfy the scaling relation Eq. Hll|) shows that the upper 
critical dimension of the system is 2. 



C. Solution for d = 

In zero dimensions, the problem may be solved for in 
a straightforward manner (see Appendix ■ The expo- 
nents are independent of since there is no diffusion. In 
this case, 



= 1, 
1 

2' 
= 1, 



(27) 
(28) 
(29) 



for all values of \i. Using Eqs. 10 and 116(1 . we obtain 



= 1, 
7=1. 



(30) 
(31) 



VI. EXPONENTS FOR ARBITRARY fj, AND d 

The question remains as to what the values of the ex- 
ponents are in arbitrary dimensions. They can be de- 
termined with the help of an assumption. We make the 
assumption that the critical exponents for a given /i are 
a monotonic function of dimension d. This assumption 
is reasonable as known exponents for most systems at 
their critical point (for example, the Ising model) have 
this property. 

Consider first the exponents when fj, = 0. Notice that 
the exponent <f> for fi = takes on the same value in the 
mean field or d = 2 (see Eq. (|24|l) as well as in d — (see 
Eq. (|29|) '). Assuming that (f>(d) should be monotonic in 
d, we obtain 



= 1, 



/! = 0. 



(32) 



Consider now the exponent 7 (as defined in Eq. dJ). It 
takes the value 1 in d = (see Eq. I31|) and in the mean 




FIG. 1: The variation of P(m) with m is shown at the critical 
point (bottom curve) and in the growing phase (top curve). 
The results are for fi — 0.25. A best fit for the curves give 
exponent values r c = 1.58 ± 0.04 and r = 1.25 ± 0.01. 



field limit (see Eq. I26|) . Hence, 7 = 1 in all dimensions 
or 



d + 2 



/j, = 0. 



(33) 



Numerical simulations of the /i = model in one di- 
mension had 4> ~ 1-01 an d ~ 1-47 jl] consistent with 
the above results. Knowing 0, the exponents r c and a 
can be solved for from Eqs. © and l|ll|) to yield 



d 2 + 6d + 4 
2(d + 2) ' 
d + 2 

~, 7: M = 

d + 4 



M = 0, 



(34) 
(35) 



Specializing to d = 1, the exponents reduce to r c = 1 1/6 
and a = 3/5. Again, these values are very close to the 
numerical values of 1.83 and 0.61 obtained in Ref. 0- 

When /i > 0, we can calculate the exponent t c and a as 
follows. Consider the exponent (3 defined by (m 2 ) ~ t@ at 
the critical point. Clearly (3 = a(3 — r c ). Note that when 
fj, > 0, P = 1 in the mean field analysis as well as in zero 
dimensions. Thus, using the argument of monotonicity 
of exponents, we obtain 



a(3 — r c ) = 1, fi > 
Solving Eqs. (|11|) and l|36|) . we obtain 
(3 - y)d + 2 



d + 2 
d + 2 

^Ti' 0< ^ <2 



< n < 2, 



(36) 

(37) 
(38) 



When [i — 1 or when d = 2, the results match the exact 
results derived in Sec. IIVI The exponent <\> for /i > is 



5 




FIG. 2: The variation of P(m) with m is shown at the critical 
point (bottom curve) and in the growing phase (top curve). 
The results are for fi — 0.50. A best fit for the curves give 
exponent values r c = 1.47 ± 0.04 and r = 1.17 ± 0.01. 



still undetermined, and there seems to be no easy way to 
calculate it. 

The analytical results are now compared with results 
from monte carlo simulations in one dimension. When 
fj, = 0, simulations were done in Ref. [g. As pointed 
out earlier in this section, the numerical values of the 
exponents are in close agreement with that obtained in 
this paper. We therefore concentrate on non-zero values 
of fi. The exponent that is determined numerically is t c 
for [i = 0.25 and = 0.5. The simulations were done 
on a one dimensional lattice of size 1000 with periodic 
boundary conditions. P(m) was obtained by averaging 
over 10 s realizations. 

In Fig. ^ the results for [i = 0.25 is shown. P(m) is 
measured for q = q c 0.380 when p — 1.0, and for the 
growing phase in which p is set to zero. The critical q c 
was fixed to be that value of q at which the distribution 
changed from an exponential to a power law. A best fit 
gives t c = 1.58 ± 0.04 and r = 1.25 ± 0.01. These should 
be compared with the analytical results t c = 1.583... 
and t = 1.25. 

Fig. |2 is as in Fig. ^ except that /x = 0.5 and 
q c w 0.448. A best fit gives r c = 1.47 ± 0.04 and 
t = 1.17 ± 0.01. These should be compared with the 
analytical results r c = 1.5 . . . and r = 1.166 . . .. 



VII. 



CONNECTION TO RELATED MODELS 



In this section, similarities between the In-Out model 
and other related models of aggregation are discussed. A 
model that closely resembles the model studied in this 
paper is the charge model with adsorption 0, 0, 0] . 
In this model, there is no longer the restriction that the 
masses have to be non-negative. Also, +1 and —1 masses 
are input at the same rate. In the limit of large time, 



P(m) for this model is a power law P(m) ~ |m| -5 / 3 for 
|m| 3> 1 in one dimension. This model could be expected 
to have the same behavior as the In-Out model at the 
critical point, since the growth velocity is zero. However, 
the exponent for the charge model is different from the 
value 11/6 obtained for the In-Out model, showing that 
the restriction of non-negative masses is relevant. We 
now ask whether it is relevant when fi > 0. 

Since the charge model was not studied earlier for 
fi > 0, we now give a short derivation of the power law 
exponent using scaling arguments. It was shown based 
on very general arguments that for the charge model, 
irrespective of the diffusion rates that 0] 



t, 



(39) 



Assuming scaling for P(m, t) as in Eq. I|12|) and the scal- 
ing relation Eq. i|ll[) . it is easy to see that P(m) ~ mT Tch , 
where 



Tch 



(3 - y)d + 2 
d+2 ' 



< (i < 2 



(40) 



But this is the same as r c obtained for the In-Out model 
(see Eq. (|37ll . Thus, the charge model and the In-Out 
model appears to have the same behavior when fi > 0. 

A reason for this could be the following. When n = 0, 
there is a chance that large positive masses get neutral- 
ized by large negative masses in the charge model. This 
process is totally absent in the In-Out model, resulting 
in the exponents being different. When fi > 0, large pos- 
itive and negative charges get immobilized and their col- 
lision becomes infrequent. Hence, one could ignore this 
process in the charge model and hence the two models 
become qualitatively similar. A pitfall of this argument 
is that it predicts that t c for /i — should be less than 
T c h, contrary to what is seen. Thus, the exact connection 
between these two models remains unclear. 

Another model which is related to the In-Out model 
is a model of coagulation with fragmentation 0, l2fl l2ll | . 
In this model, the desorption at a site is accompanied 
by adsorption at the neighboring site, thus conserving 
mass locally. In this model, there is a phase transition 
from a phase in which P(m) is exponentially distributed 
to one in which it is a power law, accompanied by a 
infinite aggregate which accommodates a finite fraction 
of the total mass. When the diffusion constant is mass- 
dependent, it was shown that pi), the exponent t c in the 
mean field limit is exactly the same as that of the In-Out 
model in the mean field limit. However, in dimensions 
lower than the upper critical dimension, the exponents 
in the model with mass conservation remains equal to 
the mean field value, unlike the In-Out model. 



VIII. 



SUMMARY AND CONCLUSIONS 



To summarize, the exponents characterizing the phase 
transition from a phase with finite mean height to one 
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in which it grows linearly with time in the In-Out model 
were calculated. The model was extended to one in which 
particles diffuse with a mass dependent rate D(m) oc 
Tji ^ . The exponents were shown to have a discontinues 
jump at \i = 0. The exponents are unrelated to previ- 
ously studied universality classes of nonequilibrium phase 
transitions. 

There are several questions that remain unanswered. 
Other models which show a wetting transition as seen 
in the In-Out model have exponents which can be ex- 
pressed in terms of exponents of absorbing phase transi- 
tions [UEi. Here, there seems to be no apparent con- 
nection to any underlying absorbing phase transition. It 
would be interesting to find connections to other models 
of nonequilibrium phase transitions. 

The calculation of exponents in this paper for arbitrary 
p relied on the assumption that the exponents are mono- 
tonic with dimension. While simulations do support the 
results that are obtained, it is important to have a more 
rigorous derivation of the exponents without making this 
assumption. Also, one would expect logarithmic correc- 
tions to the power laws in two dimensions. These have 
been ignored in this paper. A calculation of these correc- 
tions would be of interest. 

A connection to the charge model was pointed out in 
Sec. lVllI The models seem to be similar for p > 0, while 
different for p = 0. The precise connection between the 
two would be worth exploring since the charge model is 
analytically more tractable. Finally, the discontinuity of 
exponents at p — remains a puzzle. 

APPENDIX A: EXACT SOLUTION FOR p = 1 

In this appendix, we derive the exponents for p = 1 
in arbitrary dimensions. We do so by examining the two 



where C(x) = (m^mo), C M (x) = (m x mj M ), p = (m) 
and £>(x) = (m x 5 mOi0 ). 

Consider Eq. (|A6|I in the steady state when the time 
derivative is set to zero. For arbitrary values of p, the 
right hand side of Eq. (|A6|I involves three unknowns: 
C M (x), Z?(x) and p. However, when p = 1, a simplifi- 
cation occurs because 



point correlations in the system in the steady state. To 
fix notation, let x' denote one of the 2d nearest neighbors 
of the site x. Let ?y(x, x', t) be the mass transferred from 
site x to x' at time t in a time interval At. From the 
definition of the model, it follows that 

77(x,x',<) = ( m * with P rob -2^< (Al) 
1 otherwise. 

To order At, the only nonzero two point correlation in 
the noise is 

(?7 (x,x') 2 > = ^^. (A2) 

Let J(x, t) be the mass transferred due to adsorption and 
desorption from the site x at time t in an infinitesimal 
time At. Then, 

{1 with prob qAt, 

-1 + <W,o with P rob P At > ( A3 ) 
otherwise. 

To order At, the only nonzero two point correlation in 
the input / is 

(I(x,i) 2 } =(q + ps)At, (A4) 
where s = ^] , P(m) is the occupation probability. 
The mass m x (i) at lattice site x at time t evolves as 

m x (i + Ai) = m x (t)-^7/(x,x',t) + ^7 ? (x',x,t) 

x' x' 

+ J(x,i). (A5) 

To obtain the two point correlations, we multiply m x (t + 
At) by mo(t + At) and take averages over the possible 
stochastic moves. Using Eqs. (|Al|) - (IA5p . we obtain 



(A6) 



thus reducing the number of unknowns to two. Define 
F(k)=£[£(x)-p(l- s )K k - x . (A8) 

X 

Solving for F(k) from Eq. i|A6(l . we obtain 



dC*(x) 



1 d 

= -2C ll (x) + -J2 C fl (x 1 ,...,x j + k,...x d ) + 2(q-p)p + 2pD(x) 



dt " w d 

3=1 k=±l 

1 d 

+S x ,o(q + ps) + (m 2 M )(2(5 Xj0 - ^,0 • • • 

j=l k=±l 



r 



Ci(x)=p-Z>(x), 



(A7) 



(A9) 
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where 



APPENDIX B: MEAN FIELD SOLUTION 



h(k) = ^cos(fci) - (1 +p)d. 



(A10) 



To obtain p, we use the fact that the constant term in 
F(k) equals -p(l - s). Then 



dpqgjp) 
dp 2 g(p) - 



(All) 



where 

g(p) 



~2~7 



dkrl 



2^ (l+ P )d~J2ti c os(h) 



Thus, the mean density diverges as (q — q c ) , where 

q c = dp 2 g(p). (A13) 
Specializing the result to one and two dimensions, 

n 3/2 



q\° = 



„2D 



P 



2p 2 
tt(1+p)' 



K 



1 



1+p 



(A14) 



where K is the complete elliptic integral of the first kind. 

The form of the divergence of the density as q ap- 
proaches q c , namely (to) ~ q^ 1 means that 



0(2 - Te ) = 1, for p = 1, 



(A15) 



in all dimensions. To obtain one more relation between 
the exponents, we consider d< "™ ^ for large times at the 



critical point. Firstly, we need to invert Eq. (|A9I) to ob- 
tain -D(l) where 1 denotes the site (1, 0,0,...). Inverting, 
we obtain 

D(l) = fazA + (g _ pp ) [i _ d (i + p ) 5 (p)] . (A16) 

We now make the assumption that the leading time de- 
pendence to D(l) is obtained by restoring the time de- 
pendence of p. Then, when x = 0, Eq. 1|A6(1 reduces 
to 



d(r 



dt 



2p(t) {1+P)iq gc) +2qd(l +p)g(p). (A17) 



P 



If we now take the limit q — > q c before t — > oo, then we 
obtain that (to 2 ) ~ t. Thus 



a(3 



1, for /i = 1, 



(A18) 



in all dimensions. Solving for the exponents r c , a and 
from Eqs. (fTTf) . (|AT5|) and (|AT8)l . we obtain 



2d + 2 
d + 2 ' 
d + 2 



a 



(A19) 
(A20) 
(A21) 



Correspondingly 6 = (d + 2)/2 and (3 = 1. 



In this appendix, we derive the exponents for /i > 
using a mean field approximation. This approximation 
involves ignoring correlations between the masses, i.e., re- 
placing joint probability distribution functions by prod- 
uct of single point distributions. Then, the master equa- 
tion for the temporal evolution of P(m, t) is 



dP(m) 
dt 



dP(0) 
dt 



= -(m ' l + s' +p + q)P(m) + qP(m - 1) 



pP{m + 1) + 53 



P(m')P(m — m') 



= ss 1 — q + qs + pP(l), 



(Bl) 
(B2) 



where s = 53m=l P{ m ) ancl s ' = Em=i m~^P(m). The 
different terms enumerate the number of ways the mass 
at a certain site can change. Then it follows that 



d(m n 
dt 



7i — k — fj, 



fe=l 

EfDfe+c- 1 )^]^"-' 
fe=i 



) + g+(-l)> 



(B3) 



Consider first the steady state solution of Eq. (|B3|) when 
the time derivatives may be set to zero. Then, putting 
n = 1, and solving for the occupation probability s, we 
obtain 



q 

P' 



(B4) 



The mean field equations take on a simpler form for 
the cases fi = and p = 1 , and hence we solve them sep- 
arately from the arbitrary p case. Though all the mean 
field exponents for p = were derived in Ref. @ using 
the generating function method, they will be rederived 
here using a different method which will be simpler to 
generalize to the p > case. 



1. p = 



On choosing n = 2 in Eq. I|B3|) and taking the steady 
state limit, a quadratic equation for (to) is obtained 
which can be solved to yield 



p-q - \/(p - q) 2 - 4g 



(B5) 



The expression for (to) becomes invalid when expression 
under the square root sign becomes negative, thus fixing 
q c . Solving, we obtain 



9< 



= p + 2-2yi+^, 



(B6) 
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where the sign is chosen by the condition that q c = 
at p = 0. Consider now the equation corresponding to 
n = 3 in Eq. (|±33|) . Solving for (to 2 ), we obtain 



2. /Lt = 1 



(p + g)W 

p — q — 2 (to) 



(B7) 



But the denominator tends to zero as ^/q, where q = 
q c — q. Therefore, near the transition point (to 2 ) diverges 
as 



(TO 



1 

W 2 ' 



(B8) 



Consider now the equation corresponding to n = 4 in 
Eq. I|B3|I . Solving for (to 3 ), we obtain 



(to 3 ) 



3(m 2 ) 2 + 3{q + p)(m 2 )+2(q-p)(m)+q 



1 

~372> 



2(p-«-2<m» 
9^0. 



<B9) 
(BIO) 



Knowing the behavior of (to 2 ) and (to 3 ) near the tran- 
sition point, we immediately obtain 



0(3 - 7-c) 

0(3 - r c ) 
which can be solved to give 



1 

2' 

3 

2' 



5/2, 
1. 



(Bll) 
(B12) 



(B13) 
(B14) 



To calculate and r we had first taken the limit t — > oo 
followed by the limit g — * 0. In order to calculate a, we 
need to take the limits in the opposite order, namely 
q — > followed by t — > oo. We first note that in this limit 
s = q/p and (to) = (p — q)/2. Then choosing n — 3 in 
Eq. IB3I we obtain 



d(TO 3 ) _ 3(p 2 - qp 
dt ~ 2 

Thus, 

a(4 - t c ) = 1. 
Substituting for r c , we obtain 

2 



(B15) 



(B16) 



(B17) 



Thus, in the mean field limit, the scaling function takes 
on the form 



, _ . 1 / _ to \ 



At = 0. (B18) 



We start again with Eqs. i|B3|) and (|B4|I . For // 
1 . Eq. (|B3|I simplifies because (m"~ fc ~ A1 ) reduces to an 
integer moment of to. Choosing n = 2 in Eq. I|B3|I in the 

steady state, we obtain 



pq 



Tr-pq-q 

The mean mass (to) diverges when q = q c , where 

,2 



1c 



P 



l+p 



(B19) 



(B20) 



Choosing n = 3 in Eq. i|B3l) in the steady state, we obtain 
{m) 2 + (q + p){m) 



[m =p- 2 

p -pq- q 



(B21) 



Equations (|B19|I and (|B2 1|) imply that (to) ~ q 1 and 
(to 2 ) ~ q~ 3 when q — > 0. Thus, 



Solving, we obtain 



0(2 -T ) = 1, 

0(3 -t c ) = 3. 



2 
2. 



(B22) 
(B23) 



(B24) 
(B25) 



In order to calculate a, we need to set q = q c and take 
the large time limit. Then choosing n = 3 in Eq. IB3I we 
obtain 



d ( m2 ) . t 2 Qc (2-r)-l + 2j3 2 



(i+ P y 



(B26) 



If we assume that the first term is the dominant term, 
then we reach a contradiction for t (namely, r = 1). 
The other alternative is to that (to 2 ) ~ t, implying that 
a(3 — t c ) = 1 or 



(B27) 



Thus, in the mean field limit, the scaling function takes 
on the form 

P(m,Q,t)~ -^^(mg 2 ,-^), M = 1. (B28) 



3. A* > 

We will follow the same procedure as for the fi = 
and the [i = 1 cases. However, for arbitrary (i, we are 
no longer able to determine neither the critical value q c 
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nor the exponent 0. Consider the equations arising from 
choosing n = 2 and n = 3 in Eq. (|J33|> in the steady state: 



(B29) 
(B30) 



(m 1 -") = p-q-* 
(to) 



/ 2-m = 



- (p + q). 



To satisfy Eq. i|B30(l . we require that (m) diverges at the 
critical point. Substituting the scaling form, we obtain 

0(3 — A* — Tc) = 0(3 - r c ) - 20(2 - r c ) (B31) 

implying that 

r c = 2-|. (B32) 

To calculate the exponent a, all we require is that 
(to 1_ai ) is finite at the critical point and is equal to p — q. 
This follows from Eq. (|B29|) . Now, if we stay at the tran- 
sition point, we obtain that d(m 2 )/dt = 2q implying that 



t. 



Therefore a(3 — r c ) = 1 or 



2 + fj, 



(B33) 



(B34) 



Thus, in the mean field limit, the scaling function takes 
on the form 



< n < 2. 
(B35) 



APPENDIX C: SOLUTION FOR d = 

In zero dimensions, the problem becomes analytically 
tractable as diffusion no longer plays a role. Therefore, 
the exponents are independent of p. The master equation 
for the evolution of the mass distribution P(m) is 



dP{m) 

dt 
dP(0) 

dt 



= -(p+q)P(m)+pP(m+l)+qP(m-l%Cl) 
= -qP(0)+pP(l). (C2) 



The steady state solution is obtained by setting the time 
derivatives to zero. It is then straightforward to obtain 



P{m) = 



P 



q i q 

v 



TO > o. 



(C3) 



This solution is valid when q < p. For q > q c = p, there 
is no nontrivial steady state solution. The typical mass 
diverges as q approaches q c as (q — <z c ) _1 ; therefore 

0=1. (C4) 
Also, the occupation probability s = 1 when q = q c . 
Since s cannot increase beyond 1, it remains stuck at 1 
for all further values of q. When q = q — q c is positive, 
d(m)/dt = q — ps = q. This means that 9 = 1. 

The exponents r c and a may be obtained by solving the 
problem at q — q c . In this case if one were to identify to 
as the coordinate of a random walker, then the problem 
reduces to a problem of a random walker with a reflec ting 
barrier at the origin. This problem is easily solved [22] 
and in the limit of large time, 



P(m,t) 



' : expf^), t»l. (C5) 



'nqt 



4qt 



The exponents r c and a may be read off from Eq. (|C5|) 
to be 



a 



1, 
1 

2' 



(C6) 
(C7) 



in zero dimensions. 



When q > q c , the problem reduces to the problem of 
a random walker with a drift and a reflecting barrier at 
the origin. Again, this problem is easily solvable: 



„. „ 1 TO „ / TO 

P(m,t) = —- -6 



- 1 , m,t » 1. 



m(q-p)t \(q-p)t 

(C8) 

Clearly, the exponent 7 = 1 (see Eq. lfTS)l for definition). 
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